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ABSTRACT 

We obtained Herschel PACS and SPIRE images of tiie tiiermal emission of the debris disk around the A5V star/? Pic. The disk is well resolved 
in the PACS filters at 70, 100, and 160 //m. The surface brightness profiles between 70 and 160 //m show no significant asymmetries along the 
disk, and are compatible with 90% of the emission between 70 and 160/jm originating in a region closer than 200 AU to the star. Although 
only marginally resolving the debris disk, the maps obtained in the SPIRE 250 - 500pm filters provide full-disk photometry, completing the 
SED over a few octaves in wavelength that had been previously inaccessible. The small far-infrared spectral index (/3 = 0.34) indicates that the 
grain size distribution in the inner disk (<200 AU) is inconsistent with a local collisional equilibrium. The size distribution is either modified by 
non-equilibrium effects, or exhibits a wavy pattern, caused by an under-abundance of impactors which have been removed by radiation pressure. 
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1. Introduction 



The j8 Pic disk, discovered by IRAS dAumann et alJll984l) . was 
the first debris disk to be directly imaged in scattered light 
(ISmith & Terrild Il984t) . It is seen close to edge-on and ex- 
tends in the optical ou t to 95", coiTesponding to 1800AU 
dLarwood & Kalasll2001h . 

/?Pic (A5V) is one of the closest (19.44 + 0.05 pc. 



Ivan LeeuwenI 2007) and youngest d ebris disks. The estimated 



)ungi 



age (12 Mvr. 'Zuckerm an et alJl200 Jh significantly exceeds typ- 
ical timescales for th e surviv al of pristine circumstellar dust 
grains (e.g., Fedele et al . 2010), hence continuous replenishment 
of the dust, presumably through collisions of planetesimals, is 
needed. The closeness of the object e nsures that it can also be 
spatially resolved at long wavelengths: Holland et al. (1998) re- 
solved the disk at 850//m and Liseau et al. (2003) at 1200 yum. 

Optical and near-infrared observations of the inner part 
(<100AU) of the disk yield evidence of asymmetries such as 
warps and density co ntrasts, which may relate to the pres- 
ence of planetesimals (Kalas & Jewitt 1995; Pantin et al. 1997; 
Mouillet et al. 1997; Heap et al. 2000 ; Telesco et al. 2005). 



Lagrange et alJ (l2009h imaged a possible companion at a pro- 



jected distance of 8 AU from the star. 

Images of jSPic in scattered stellar light directly detect small 
grains and indirectly larger grains that produce the smaller ones 
through collisions. The grain-size distribution can be quantita- 
tively constrained from the spectral energy distribution (SED) of 
the disk, in the infrared and (sub)mm domai ns. The spectral in- 
dex of the SED at the longest wavelengths dLiseau et al. 2003; 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



iNilsson et al.ll2009l) is inferred to be fairly low, which according 
to modeling (fcraine 2006; Natt a et al-IIJOOTl; iRicci et alJl2010l) 
can be interpreted as a deficit of small grains. 

In this paper, we present far-infrared imaging of the /? Pic 
debris disk in six Herschel photometric bands between 70 and 
500 yum. These bands cover for the first time the long-wavelength 
side of the peak in the thermal emission of the disk, and the large 
aperture of the telescope enables us to resolve the disk at far-IR 
wavelengths for the first time. With these data, we measure the 
surface brightness profiles of the disk and readdress the issue of 
the grain-size distribution in the inner 200 AU. 

2. Observations and data reduction 

We obtained maps o f ISPic with the PA CS and SPIRE instru- 
ments of Herschel (Pilbratt et al. 2010). The in-orbit perfor- 
mance, scientifi c capabilities, c a librati on methods, and accuracy 
are oudi ned bvlPoglitsch et al.1 (1201 Oh for the PACS ins trument 
and by Griffin et all (1201 Oh and lSwinvard et all ([2010^ for the 
SPIRE instrument. The observations were carried out during the 
science demonstration phase as proposed in the 'Stellar Disk 
Evolution' guaranteed time proposal (PI G. Olofsson). Table lA. fl 
gives a summary of the observations. The deep PACS observa- 
tion at 70/im and 160/im is a standard PACS photometer scan 
map, split into a scan and cross-scan on the sky. The sky scan 
speed was 10" sec"'. The homogeneously covered area of the 
deep map is 2.5' x 2.5'. The observation at 100 //m is much shal- 
lower, with a single scan direction at a rate of 20" sec"', homo- 
geneously covering an area of 2' x 2'. The PACS beams at 70, 
100, and 160A(m are 5.6, 6.8, and 11.3"FWHM. In the SPIRE 
observation, the three bands are observed simultaneously in a 
standard scan map. The map coverage is 8' x 8'. The SPIRE 
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Table 2. Overview of measured quantities: position angle PA, 
northeast (NE) and southwest (SW) extent (signal reaching the 
1 cr noise), map noise level, beam FWHM, and the flux density 
integrated over a 60" aperture. 



FWHM beam sizes in the 250, 350, and 500 //m channels are 
18.1, 25.2, and 36.9" respectively. 

The data processing is described in Appendix. The absolute 
flux calibration accuracy of the resulting PACS maps is better 
than 10% at 70 and 100//m, and 20% at 160//m ('Poglitsc h et"an 
[2010). The flu x calibration accuracy of the SPIRE maps is bet- 
ter than 15%. dSwinvard et al.ll2010l) . The Icr noise levels of the 
maps are listed in Table |2l 



3. Analysis 

In Fig. [1] we show the maps obtained in the three PACS fil- 
ters (70, 100, and 160//m) and flie three SPIRE filters (250, 
350, and 500 yum). We also compare the point spread functions 
(PSFs) measured on the asteroid Vesta using the same satellite 
scan speed, processed as the y6 Pic maps and rotated to align with 
the telescope pupil orientation on the sky during the fi Pic obser- 
vations as listed in online Table lA.fl 

These images show a clearly resolved disk from 70-160//m. 
Each map was fitted using a 2D Gaussian function. Within the 
2" Herschel pointing accuracy, the Gaussian center matches 
the star's optical position. The fitted position angles, listed in 
Table |2l agree wi th the optical d isk position angle of 30°8 re- 
ported by iKalas & JewittI ([1995). Cross-sections orthogonal to 
the disk position angle in the NW to SE direction show no signif- 
icant broadening compared to the PSF. The disk is not resolved 
in the vertical direction. The feature towards the NW, visible 
in the 70-160/vm images, is produced by the three-lobed PACS 
PSR 

In Fig.|3j we present the surface brightness profiles along the 
disk position angle. We compare them with the cross-sections 
aligned in the same direction through the PSFs. At 250 and 
350 yum, the disk is marginally resolved. At 500 yum, the fi Pic 
profile shows no significant departure from the PSF profile, 
with the exception of a cold blob in the southwest. As can be 
seen in Fig.[Tl the location of this feature in the 250 - 500 yum 
maps coincides with the fl ux p eaks seen at 850 and 870yum 
bv Hofland et al. (199^ and lNilsson et al. (2009) . respectively. 
However, the lOOarcmin^ region around y8 Pic (depicted in on- 
line Fig IA.2I 1 shows more than 50 background sources compa- 
rable to this feature in the 250yum map. The feature is therefore 
probably a background source. 

Other asymmetries between the northeast and southwest pro- 
file are within the errors induced by the asymmetry of the PSF. 
No sharp disk edge is seen; in all filters, the surface brightness 
declines gradually to the Icr detection limit of the maps. Table|2] 
lists the extent of the detected emission region in the NE-SW 
direction. 
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Fig. 4. Normalised surface brightness profiles along the disk in 
NE direction in the three PACS filters. The profiles were con- 
volved with a Gaussian to match the spatial resolution of the 
160 jum image. The same convolved profiles are shown for the 
PSF maps. 



The comparison of the surface brightness profiles in the three 
PACS filters in Fig. |4] shows the same brightness profile along 
the 30.8° position angle in NE direction. The 70 and 100//m 
profiles were convolved with a Gaussian to match the spatial 
resolution at 160 jum. The same convolution was applied to the 
70 and 100 yum PSF profiles. The shape of these convolved PSF 
profiles defers significantly from that of the 160//m PSF profile. 
The wiggles in the 160yum profile differ up to a factor of 3 from 
the convolved 70 and 100 //m PSF profiles. Within these uncer- 
tainties, there is no evidence of a wavelength dependent surface 
brightness. This indicates that the grains producing the emission 
at 70, 100, and 160yL/m are confined to the same locus in the disk. 
At 70jum, the broadening of the profile with respect to the PSF 
indicates that 90% of the emission originates in a region within 
11" or 200 AU of the star. 



4. The far-infrared SED and grain size 

We integrated the surface brightness maps over a 60" radius cir- 
cular aperture. Background subtraction was based on a rectan- 
gular region, selected close enough to the object to be within the 
map region with the same coverage as the center of the map. 
For the background outlier rejection, the DAOphot algorithm in 
the HIPE aperture photometry task was used. The aperture pho- 
tometry obtained provides a good measure of the flux density of 
the integrated disk. The contribution of the stellar photosphere 
at these wavelengths is negligible. The error is dominated by the 
present uncertainties in the absolute flux calibration of both in- 
struments. The full disk flux densities are listed in Table|2l 

Figure |5] shows the new PACS and SPIRE photometry, and 
selected infrared and (sub-)mm flux densities from the litera- 
ture. Because the disk is optically thin at these wavelengths, 
the wavelength dependence of the emission directly probes the 
dust grains, and, in particular, their size distribution. We over- 
plot two modified Rayleigh-Jeans laws (Fy oc v^^^^*), normalized 
to the 160yum datum. The spectral index/? indicates the mean 
dust opacity k oc \fi. An index /3 - corresponds to a black 
body with a k independent of wavelength A, indicating grains 
that are much larger than A/In. Interstellar grains, which have a 
size distribution f(a)oca"* with q - 3.5 and an upper siz e Umit 
of a^ax ~ 0.3 yum, are characterized by y6 = 1.8 + 0.2 (IDraind 
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Fig. 1. Surface brightness maps of the /? Pic debris disk at 70, 100, 160, 250, 350, and 500 //m. The PACS PSFs, rotated to match 
the position angle of the telescope at the time of the ySPic observations are depicted in the upper right corner of the images. The 
SPIRE PSFs are depicted in Fig. lA.6l All images are scaled linearly, contour lines are in steps of 10% of the peak flux. The surface 
brightness unit is Jy arcsec"^. Th e white circle shows the b eam FWHM. The p osition of the flux peaks observed at 850, 870, and 
1200 A/m bv lHolland et~an (Il998l) . iNilsson et al] (l2009h . and lLiseau et"an (l2003h are indicated with H, N, and L. 



l2006l) . In protoplanetary disks, y6-values fr om 1.5 down to 
are found, depending on the disk geometry dAcke et al.ll2004l) . 
An eiTor-weighted least squares fit of a Rayleigh-Jeans law 
to the ySPic photo metry at waveleng ths beyond 160/im yields 
y6 = 0.34 + 0.07. Nils son et al.l (12009) obtained B = 0.67 from 
a yS-corrected black-body fit to the full disk SED, including mid- 
infrared photometry. The difference between both results should 
not be over-interpreted since both approaches are sensitive in 
different ways to simplifying assumptions about the temperature 
and size distribution within the disk . In any ca s e, bot h results 
consistently show a value below 0.7. Ricci et all (1201 Ol) demon- 
strate that such a low value cannot be explained with a.q - 3.5 
power law. This is a surprise insofar as the latter value is the ex- 
pected result for a population of bodies in a standard steady-state 
colUsional cascade dPohnanvil 19691) . 



The small value of (3 can be interpreted in a number of 
ways. The grain size distribution can exhibit a wavy pattern, 
caused by the absence of impactors small enough to be effi- 
ciently blown out of the disk by radiation pressure. This causes 
an over-abundance of the grains that are just bound, which means 
there are more impactors for the next larger size population. The 
reduction of this population caus es an over-abundan ce of a fol- 
lowing size population and so on dKrivov et al.ll2006h . The wavy 
size distributi on can lead to small value s of yS when measured 
in the FIR (Theba ult & Augereadl2007h . If the wavy structure 
were as strong as found in this paper for normal and weak ma- 
terial properties, it would be consistent with the small fi value 
we have measured. However, the strength and phase of the wavy 
pattern in the size distribution depend on both the grain structure 
and the eccentricity of the dust orbits in the disk. 



The grain size distribution in ySPic must be flatter than the 
q = 2).5 power law, meaning that the fraction of small par- 
ticles must be lower. Radiation pressure can push the smallest 
grains (with Fiad//^grav > 0.5) onto hyperbolic orbits, hence re- 
duce the time these particles spend in the inner part of the disk, 
which can decrease their v olume density by two orders of mag- 
nitude (iKrivov et al.ll200C)l) . The disk cannot be fully cleared of 
small particles, since it has been seen in scattered light out to 
1800 AU. The scattering grains are probably the (sub-)/im grains 
that are blown out of the inner disk, where the collisions take 
place. However, this effect only reduces the densities for grains 
of size below a few micrometers, and even fully removing these 
grains would not change (3 to the observed value. 



Alternative explanations of the small value of/? cannot be ex- 
cluded. There are indications that the grains produced in the deep 
impact experiment followed a flatter power law with ^ ^ 3.1 
(Jordaet al. 2007). Laboratory experiments illustrate that frag- 
ments produced in collisions of porous aggregat es can follow 
much flatter slopes (^ = 1 .2, iGiittler et al. 1 120091) . demonstrat- 
ing that the porosity of the colliding grains should not be disre- 
garded. 



Additional dynamical models should be developed to quan- 
tify the possible contribution of these effects to the small /? ob- 
served in yS Pic. 
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Fig. 3. Surface brightness profiles along the disk in NE-SW direction, following the 30.8° position angle. The black horizontal hne 
shows the Icr detection limit in the yS Pic maps. The yS Pic profiles are shown with the Icr errors. 
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Fig. 5. The infrared to mm SED of ySPic. The PACS (70 - 
160yum), and the SPIRE (250 - 500 yum) fluxes were inte- 
grated over a 60" radius aperture. IRAS flux densities are 
from the IRAS Poin t Sour ce Catalog. The 850yum SCUBA 
datumjlolland etlD [Illll) and the 1200 yum SIMBA datum 
(iLiseau et al.i i2003] ) are integrated over a 40" radius aperture. 
Overplotted is a Rayleigh-Jeans extrapolation of the 160jum flux 
density with a spectral index y6=0 andy6=2, and the best fit to the 
160-1200 yum data (J3=Q.34). The stellar photosphere is a Kurucz 
model for reif=9000K;log(g)=3.9 scaled to die 2MASS pho- 
tometry Kj=3.52. 



5. Conclusions 

We have presented images of the /3 Pic debris disk in six photo- 
metric bands between 70 and 500 /vm using the PACS and SPIRE 
instruments. We resolve the disk at 70, 100, 160, and 250yum. 
The images at 70-160yum show no evidence of asymmetries in 
the far-infrared surface brightness along the disk of ySPic. The 
observed profiles are compatible with 90% of the emission orig- 



inating in a region within a radius of 200 AU from the star. The 
disk-integrated photometry in the six Herschel filters provides a 
far infrared SED with small spectral index/? ^ 0.34, which is 
indicative of a grain size distribution that is inconsistent with a 
local collisional equilibrium. The size distribution is modified by 
either non-equilibrium effects, or exhibits a wavy pattern, caused 
by the under-abundance of impactors that are small enough to be 
removed by radiation pressure. 
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Table A.l. Observation log 





observation 


date 


Pos Angle 


duration 


filters 


SPIRE 


1342187327 


2009-11-30 


154.96° 


3336 s 


250,350,500 


PACS 


1342185457 


2009-10-07 


106.54° 


866 s 


100,160 


PACS 


1342186613 


2009-11-01 


130.98° 


5506 s 


70, 160 


PACS 


1342186612 


2009-11-01 


130.98° 


5506 s 


70, 160 



Appendix A: Data reduction 

The PACS data were processed in the Herschel interactive anal- 
ysis environment HIPE (v3.0), applying the standard pipeline 
steps. The flux conversion was done using version 5 of the 
response calibration. Signal glitches due to cosmic ray im- 
pacts were masked out in two steps. First the PACS photMMT- 
Deglitching task in HIPE was applied on the detector time- 
line. Then a first coarse map was projected, which is then 
used as a reference for the second level deglitching HIPE task 
IlndLevelDeglitch. In the detector time series we masked the 
region around the source prior to applying a high-pass filter to 
remove the low frequency drifts. The scale of the high pass filter 
was taken to be half the length of an individual scan leg on the 
sky, i.e. 3.7'. The detector time series signals were then summed 
up into a map using the PACS photProject task. The pixel scale 
for the 70 and 100 //m maps was set to 1", while the scale for the 
160yum map was 2". For the deep map in the 70 and 160jum fil- 
ter we combined the two detector time series and projected these 
together. 

The SPIRE data were also reduced using HIPE and maps 
were obtained via the default naiveMapper task. The SPIRE ob- 
servation consists of several repetitions of a map observation of 
the same area. As a result it was possible to project the data with 
a pixel size of 4, 6, and 9" while still maintaining complete sam- 
pling across the source. 
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Fig.A.6. The 250, 350 and 500 /im SPIRE PSFs, rotated to match the position angle at the time of the yS Pic observations. The PSF 
images are scaled linearly, contour lines are in steps of 10% of the peak flux. The white circle shows the beam FWHM. 




Fig. A.2. The 250 jum SPIRE map around the fi Pic disk. The lOxlO'region delimited by the white square shows more than 50 
background sources comparable to the cold blob seen in the southwest of the disk. 



